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Abstract 

We study multiple defaults where the global market information is modelled as 
progressive enlargement of nitrations. We shall provide a general pricing formula by 
establishing a relationship between the enlarged filtration and the reference default- 
free filtration in the random measure framework. On each default scenario, the formula 
can be interpreted as a Radon-Nikodym derivative of random measures. The contagion 
risks arc studied in the multi-defaults setting where we consider the optimal investment 
problem in a contagion risk model and show that the optimization can be effectuated 
in a recursive manner with respect to the default-free filtration. 



1 Introduction 

The contagion credit risk analysis with multiple default events is an important issue for 
evaluating the credit derivatives and for the risk management facing the financial crisis. 
Compared to the single credit name studies, there are several difficulties in the multi- 
defaults context. Generally speaking, the global market information containing all defaults 
information is modelled as a recursive enlargement of nitrations of all default times with re- 
spect to a default-free reference filtration. To obtain the value process of a credit-sensitive 
claim, one needs to consider the conditional expectation of its payoff function with respect 
to the global market filtration. The mathematical formulation and computations are in 
general complicated considering all possible default scenarios and the enlarged filtration. 
Furthermore, the modelling of correlation structures of default times in a dynamic manner 
is a challenging subject, in particular, when it concerns how to take into consideration the 
impact of one default event on the remaining names. 

In the literature, there are mainly two approaches - bottom up and top down - to 
model multiple default events. In the first approach, one is interested in the probability 
distributions of each individual default and in their correlations, often using copula func- 
tions. The second approach concentrates directly on the cumulative losses distributions, 
which allows to reduce the complexity of the problem. However, the correlation structure 
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between the default times is not straightforward in the top-down models. Recently, a new 
approach has been proposed to study the successive default events (see [5], also [3]), which 
provides an intermediary point of view between the above two approaches. The main ideas 
of [5] are two-folded. On one hand, the default scenarios are largely reduced and we can 
thoroughly analyze the impact of each default event on the following ones; on the other 
hand, the computations are decomposed on each default scenario and hence concern only 
the default-free reference filtration. One key hypothesis is that the family of default times 
admits a density with respect to the reference filtration. The density hypothesis is a stan- 
dard one in the enlargement of nitrations (see for example [7] for the initial enlargement 
of filtration). In the credit risk analysis, the density hypothesis has been adopted in [3] 
for analyzing what goes on after a default event, it has also been proved to be useful in 
the recursive bef ore-default and after-default extensions. 

Inspired by the second idea mentioned previously in [5], we study the non-ordered 
multiple defaults by establishing a formal relationship between the global market filtra- 
tion and the default-free reference filtration. We shall adopt the framework of random 
measures, which will provide us convenient and concise notations. One can find a detailed 
introduction to random measures in [8l Chap II] and in the monograph [2] . Similar notions 
have also been used in the filtering problems (see e.g. [12] and [IS]). The main advantage 
of introducing such a general framework is that we can treat the multiple defaults case in 
a coherent way as in the single name case. Another consequence is that we can remove 
the density hypothesis. In the case where explicit results are needed and where the den- 
sity hypothesis is assumed, we recover a result in [5], to which we refer for more detailed 
discussions. 

As applications, we are interested in the pricing with multiple defaults and in the 
contagion risks. The important idea in both cases, as mentioned above, is to find a suitable 
decomposition of the problem on each default scenario, so that the analysis will only 
concern the default-free filtration. We shall present a general pricing formula, which gives 
the value process with respect to the global market information. On each default scenario, 
the formula can be interpreted as a Radon-Nikodym derivative of random measures. This 
result can be applied to credit portfolio derivatives and also to contingent claims subjected 
to contagion default risks. 

The contagion credit risk during the financial crisis is an important subject which needs 
to be taken into consideration. Notably, one default event can have significant impact on 
remaining firms on the market and may become the potential cause of other defaults, 
called the contagious defaults. We shall present a contagion risk model to describe this 
phenomenon, where each asset is influenced by the default risks of an underlying portfolio 
and has a jump on its value at every default time. We consider furthermore an investment 
portfolio containing such assets and study the optimal investment strategy. We show 
that the global optimization problem is equivalent to a family of recursive optimization 
problems with respect to the default-free filtration. 

The paper is organized as follows. We present in Section 2 the mathematical framework 
of random measures. A general pricing result concerning multiple defaults is deduced using 
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random measures in Section 3, and is applied to credit portfolio derivatives. The Section 4 
is devoted to the analysis on contagion risks. We firstly present a multi-defaults contagion 
model, and then study the optimal investment problems in the presence of contagion risks. 

2 Random measure framework 
2.1 Preliminaries 

Let (fi, G, P) be a probability space equipped with a reference filtration F = (Ttjt^o 
satisfying the usual conditions. Let r = (n, • • • ,r n ) be a family of random times taking 
values in M", considered as K^ 1 '"'™ , representing the family of default times. Denote by 
the index set {1, ••■ , n}. We suppose that T{ (i G B) are strictly positive and finite, 
and that r, ^ tj, a.s. for i j (j G O). Let B* = (T>l) t >o be the smallest right-continuous 
filtration such that Tj is a ©'-stopping time. More precisely, T>\ := C\ s >o a ( T i A (i + £ ))- 
Let G = (Gt)t>Q be the progressive enlargement of F by the default nitrations, namely, 
G = F VD 1 V • • • VB'". 

For any I C 9, let 77 = (rj)j g /, which is a random variable valued in For i G E + , 
the notation denotes the event 



The events (AI)j c q describe all default scenarios at time t. Note that O is the disjoint 
union of (Aj.)i C Q. 

The following lemma is an extension of a classical result on progressive enlargement of 
nitrations with one default name. 

Lemma 2.1 For t G R+, any Gt -measurable random variable Yt can be written in the 
decomposed form 



where Y/(-) is a Tt ® B(M. I + ) -measurable function on x B(R^_) being the Borel 
a-algebra. 

Proof. By definition, for any t G and any integer m > 0, the random variable Yt is 
Tt V ct(t A (t + l/m))-measurable. Hence there exists an ^ ® £>(M™ )-measurable function 
F m such that F< (u>) = F m (cj, r A (t + 1/m)). For / C 0, 




(1) 




ice 
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So, for fixed u>, one has 

l A i(uj)Y t (uj) = \ A l{0j)F m (uj,T I {L0), (*+-)/=) 
when m is large enough. Let 

Y t \si) := hm sup i^u;,^) 

m— >oo 

where = (x[ m \ • • • , Xn ) is defined as x\ m ^ := Sj if i € I and x[ := t + 1/m if i ^ I. 
Then one has H A iY t = \ A iY t I {r I ). □ 

The following variant of Lemma |2. II will be useful further on. 

Lemma 2.2 Any Qt- -measurable random variable Yt can be written as 

Y t = J2*A{_Y t I (Ti), 



ice 



where 



4-:=(n^<*})n(n{^>*})> 



and Y/(-) is Tt~ ® B(M.+) -measurable. 

Proof. For any t > 0, Qt~ = [_J £>0 e V a(r A (t — e)). So there exists an integer 
m > such that Yt is <8> cr(r A (t — l/m))-measurable. Hence 1^ can be written as 
F m (uj,T A (t — 1/m)) where F m is an Tt~ <8> B(Wi. )-measurable function. Then we can 
complete the proof by a similar argument as for Lemma 12.11 □ 

Remark 2.3 In the Lemmas 12.11 and 12.21 if the random variable Yt is positive (resp. 
bounded), then Y/(-) can be chosen to be positive (resp. bounded). 

2.2 Random measures 

Definition 2.4 Let ^ T be the measure on (0 xl" , J-oo®B{W\. )) such that for any positive 
and (g> B(Wj. )-measurable function /1qo(')j 

(2) j h 00 {s)ii T {du,ds)=¥,[h 0Q {T)} 1 

where JF^ = |J t >o ^ and s = (si, ■ • ■ , s n ). 

The measure /i x can be considered as a transition kernel from (£l,Foc) to (R^ , B(R!Ji)) 
whose marginal on Q coincides with P. It can also be considered as the conditional law of 
r on Toq. We give below an example of [i r using the copula model in 
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Example 2.5 (Schonbucher and Schubert) For any % G 0, define the default time by 
T{ = ini{t : A\ > U{\ where A* is an continuous increasing F-adapted process and E7j is 
an exponential distributed random variable independent of J-^. The conditional survival 
probability is u\ = P(rj > t^oo) = exp(— A|). Note that H-hypothesis holds in this model, 
that is, P(Tj > t\J-t) = P(ri > i^oo)- The construction of joint survival distribution in 
|14j is by introducing a copula function C : W] — > M + such that 

P(T > S^oo) =P(n > Sl,-- - ,T n > S„|^oo) =C«,-- - ,<J. 

Then for any positive and ^"oo <g> B(R")-measurable function /ioo(-)> one has 
E[/ioo(T)] = / ^ 00 (a) M r (dw,da)=Ef / M*)(-l)Xi • • • d fl „C«, 

J L JR™ 



where d si ■ • ■ d Sn C(u\ , • ■ ■ , u™ n ) is an n-dimensional Lebesgue-Stieltjes measure associated 
to <?«,••• ,<J. 

Classically the random measure is a straightforward extension of the notions of in- 
creasing processes and their compensators, see [H Chap II] for details. Here the random 
measure /i T is useful to define auxiliary measures on suitable c-fields. For t > 0, let [ij 
be the restriction of [i T on Tt ® B(W±). It represents the conditional law of r on For 
this reason, we also write fjj as E[^" r |^ r t]- If ht(s) is a positive T± <S) B(Wf )-measurable 
function, then 

E[/i t (r)] = J h t (s)fi T (duj,ds) = J h t (s) [i[ \duj , ds) . 

For / C {1, • • • , n}, let \x\ be the measure on (fi x R+,Ft ® which is a partial 

marginal measure of //J" such that for any positive and Tt <8> £>(M+)-measurable function 
ht(si), si = (si)i e /, one has 



(3) / ht{si)nl(du),dsi) = / / h t (si)fj,J(duj,ds). 
This relation can also be written as 

(4) [i\ {duo, dsi) = / [i[ (doj , ds) . 

J]t,oo[ IC 

For T > i, J C {1, ■ • • , n} and Yt(0 which is positive and ^"t ® B(R!J. )-measurable, we 
define /i^ T,/ as the weighted marginal measure on (Q x R^_,.Ft ® B(R^_)) such that 

(5) / h t (si)^ T)I (duj,dsi) = I I h t (si)Y T (s)^' r (duj,ds). 
Similarly, we write /^ T,/ as 

(6) n] Tj {duj,dsi)= [ E[Y T (s)n T \F t ](duj,ds), 

J]t,oo\ IC 
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where E\YT{s)^i r \J r t\ denotes the restriction of the measure Yt{s)^ t on J-t®B(WL). Note 
that one has E[Y T (s)fj, T \J 7 t ] = E\Y T {s)n T \F t \. 

We shall use the Radon-Nikodym derivative of random measures to interpret diverse 
conditional expectations. 

Proposition 2.6 Let T > t > 0. For any positive and Tt <8> B(W±) -measurable function 
Yt(-) on x R™ , one has 

L n E[Y T (s)^\T t }(dio,ds) 

(7) E[F T (r)|^] = — — 

J R „ fi T {duj,ds) 

where the Radon-Nikodym derivative is taken on J~t- 

Remark 2.7 Note that J Rn fi T (dLV,ds) = F(du>), the above result can be written as 
E[Y T (r)\T t ]F(dio) = [ E[Y T ( Sl )^\F t }(du;,ds). 

Proof. Let ht be a positive .Fj-measurable random variable, then 

h t (io)K[Y T (T)\T t }F(dio)=E[h t E[Y T (T)\T t ]} = E[h t Y T (r)] = f h t (io)Y T (s)^(du,ds). 
Hence the equality holds. □ 

Remark 2.8 In particular, the conditional expectation E[lr|-7"t] where It is a positive 
(/^-measurable random variable can be written in a decomposed form. In fact, by lemma 
12.11 one has 

Y T = E ^A'YKti) = E a [0,TFx]T,oor(T)^(^), 

ice ice 
where Y^(-) is positive and Tt ® £?(M+)-measurable. Hence Proposition 12.61 gives 

E[Y T (T)\F t ]F(cLu) = y] [ E[Y T ( Sl )^ t ]{dw,ds). 



ice 



3 Pricing with multiple defaults 

For the purpose of pricing, let us consider a contingent claim sensitive to multiple defaults 
with the payoff function Yt(t) where Yt(s) is a positive and Tt ® £>(M™ )-measurable 
function on x M", T being the maturity. Since r = (ti,--- ,r n ) represents a family 
of default times, Yt(t) can describe a large class of financial products such as a basket 
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credit derivative, or a single-name contingent claim subjected to the default risks of mul- 
tiple counterparties, or a basket European option with contagion risks etc. The price of 
this product is computed as the expectation E[lj , (r)] under some risk-neutral probabil- 
ity measure. The dynamic price process given all market information at time t < T is 
the conditional expectation E[Yr(r)|^]. In this section, we shall present the evaluation 
formulas using the random measures. 



3.1 General pricing formula 



We suppose in this section that the (conditional) expectations are taken under some risk- 
neutral probability 

Theorem 3.1 Let T > t > 0. For any positive and Tt <8> £>(M") -measurable function 
Yt(-) on $7 x , the measure pj 7 ^' 1 is absolutely continuous with respect to \x\. Moreover, 
the following equality holds 



(8) 



Y T ,I 



nY T (T)\g t ] = ^2l A{ ^^(u,rj). 
ice afI t 



Using the notations (jH) and Q, the above equality can also be written as 
(9) 



e[y t (t) \g t ] = > i A i 1 1 „ — — — 



U=Tl 



where the Radon-Nikodym derivative is taken on the a-algebra 

Proof. By definition (|3|), one has that, for any M G Tt <8> /uf (M) = if and only if 

\m(uj, s/)l]i. +00 [-r c (sjc) = 0, /i T -a.e. Hence this implies pJ T,I (M) = 0. 
On the set Aj., any ^-measurable test random variable can be written in the form Zt(ji), 
where Z t (si) is positive and Tt <8> £?(i?+)-measurable. To prove (jHJ), it suffices to establish 

r du^' 1 
E[l M Z t (r I )Y T (r))=E[l Al Z t (T I )-^(u;,T I ) 



(10) 
One has 



E 



da YrJ 



Y T ,I 



%t]^( g /)%+oop(sj0^( g /) {^^^^{duj^ds) 



^[o,ty{ s i) z t{si) d ^ j (uj, si)f4(du, ds 7 ) 
afj, t 

/ 1 [o,t]'( s i) 1 ]t,+oo[i c ( s i' : ) z t(si)YT(s)p' r (duj,ds), 



Y T ,I 
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where the first equality comes from the definition of fi T , the second one comes from 0, 
the third one results from the definition of Radon-Nikodym derivative, and the last one 
comes from ([5]). Again by the definition of the last formula equals the left side of (jlOl) . 
□ 



Remark 3.2 The following form of Theorem 13.11 will be useful. Let Yt be a positive 
C?T-measurable random variable, which is written as 



T 

ice 



where Y^(-) is Tt <8> -8(R+)-measurable. Theorem 13.11 implies that 

I]t,oo[ JC E [ 1 [o,TKx]T,ooK c ( s ) y T( s /)Mf \Ft](dv, da) 



T 

ice 



E E i« 



icejce 



I]tM JC Vl(duj,ds) 



x ^ iir, oo[ ^x]t,T]A./ ]E [^r( s j)Mrl-^](^ ; ds ) 

L rJC u7(dLJ,ds) 
Jce idj J]t,oo[ J rtK ' ; 

where the last equality comes from an interchange of summation. 

Inspired by [5], we consider the case where the density hypothesis holds. Let v T be 
the marginal measure of [i r on B(W±_), that is, 

u T (U) = n T (n x U), VUeB(Rl). 

Note that v T is actually the law of r. 

Assumption 3.3 We say that r = (ti,--- ,r n ) satisfies the density hypothesis if the 
measure ^ T is absolutely continuous with respect to P ® v T . We denote by at(-) the 
density of /i T with respect to P <g> v T on (0 x R™ , T t ® 5(R+)), where f € R+. 

Under the above density hypothesis, one has, for any positive Borel function on R™ , 

E[f(r)\T t ] = [ f(s)a t (s)u T (ds). 
This relationship can also be written as 

(12) nl (du, ds) = a t (s)P(duj) ® v T (ds). 
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Corollary 3.4 We keep the notation of Theorem \3.1\ and assume in addition the density 
hypothesis \3.3l Then 

(13) E [Y T T )[ft = ^ l A ! J -^±- —7 n 



s/=rj 



Proof. By the density hypothesis (|12p and by (j3|), for any I C O, one has 



fil(duj,dsi) = / nJ(dco,ds) = / at(s)P(cL>) ® v T (ds). 

J]t,OD[ IC J]t,O0[ IC 



Similarly, 



Ht T ' I (du,ds I ) = [ E[Y T ( Sl )fj,Z\Ft](du,ds) = [ E[y T (s/)a T (s)|^]P((M^ x (ds) 

J]t,oo[ IC J]t,oo[ IC 



So 

d^ 1 I ]t ,oor nYT(s)a T (s)\^ T (ds) 
M Jjt.oop a t {s)^{ds) 

In fact, for any positive Tt <8> i3(M+)-measurable function /it(-); one has 

w J ]t , +ODr nY T (s)a T (s)\T t ]^(ds) 
ht{si) = - f . . . n t {(L},dsi) 



r L ,< IC ¥,[Y T {s)a T {s)\Tt]v T {ds) 

L, oo[/c E[y T ( s ) Q!T ( s )|^] i /-(d s ) 



/ ^ 1 \ tu \ / at(s)z/ x ((is) 



I P(du) [ h t ( Sl )K[Y T (s)a T (s)\^ T (ds) 

Jn Jw x]t,oo[ IC 

h t (s I )Y T (s)a T (s)F(duj) ® z/^ds) 

nxK^.x]t,oo[ /c 

h t (si)fj, t T ' (<ko,dsi), 



where the first and the last equality come from the density hypothesis, the second and the 
fourth ones come from Fubini's theorem. Thus the equality follows from Theorem 
O □ 



3.2 Pricing of credit portfolio derivatives 

We now apply the previous pricing formulas to the two important types of credit portfolio 
derivatives: the A: th -to-default swaps and the CDOs. 
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3.2.1 Basket default swaps 



A fc th -to-default swap provides to its buyer the protection against the k th default of the 
underlying portfolio. Let (r^)ke& be the ordered set of r = (Tj)i G ©, that is, vm < • • • < 
77 n ) . The protection buyer pays a regular premium until the occurrence of the fc th default 
time m) or until the maturity T if there are less than k defaults before T. In return, the 
protection seller pays the loss 1 — R(ty where is the recovery rate if < T, and 
pays zero otherwise. So the key term for evaluating such a product is the indicator default 
process l{T (fe) <T} with respect to the market filtration Q t . 



Proposition 3.5 For any t <T, 



(14) 



Efl 



{ T (k)> 



|J|<fe lDJ,\I\<k 



}t,oo[ J 



fif(du), ds) 



SJ = TJ 



Proof. Observe that l{ T(fe)>T }. = Yl\i\<k ^A ! T • By Theorem 13.11 one obtains 



Jd 



SJ=TJ 



By taking the sum over I such that |/| < k and by interchanging the summations, one 
gets (HU). □ 



Remark 3.6 Among the basket default swaps, the first-to- default swap is the most im- 
portant one. In this case, k = 1. Proposition 13.51 leads to 

Jir,oo[*^ T (^,<k) E[l { >r} |^] 

The last equality is a well-known result (e.g. [1]). 



3.2.2 CDO tranches 

A CDO is a structured credit derivative based on a large pool of underlying assets and 
containing several tranches. For the pricing of a CDO tranche, the term of interest is the 
cumulative loss of the portfolio It = J27=i Ri^-n<t, R% being the recovery rate of Tj. A 
tranche of CDO is specified by an interval corresponding to the partial loss of the portfolio. 
The two threshold values, a upper value a and a lower one b, defines a tranche of CDO and 
the loss on the tranche is given as a call spread written on the loss process It with strike 
values a and b. Therefore to obtain the dynamics of the CDO prices, we shall consider 
E[(l T -a)+\g t }. 

On the market, it is a standard hypothesis to suppose that the recovery rate for each 
underlying name is constant (equal to 40% in practice). We make this hypothesis below 
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and discuss further on in Remark 14.101 the case where the recovery rates Ri are random 
variables. 

The following result allows us to deduce the CDO prices using the k th -to default swaps. 



Proposition 3.7 Assume that Ri = R is constant for all i S 6, then 
(15) (l T -a)+ = R £ min(A; - -|, 1) 1 { 



k>a/R 



Moreover, 



k>a/R L |J|<fc lDJ,\I\<k J]t,oo[ J ^ V ' V 

Proof. Notice that for any m 6 0, the following equality holds 

(m-—)+= E min(fc- -|,1). 



S.J = TJ- 



R'^ ^ R 

keQ,ji<k<m 

Hence 

n 



*=i fcee,f<fc 
Since {k < Y17=i ^{n<T}} = { T (fc) < we obtain (fl~5|) and 



E[(lr - a)+|&] = £ min ( fc " l)E[l { T {fc) <T}[&]- 
By Proposition 13.51 the result holds. □ 



4 The contagion risk model 



In this section, we are interested in the contagion risks with multiple defaults. One ob- 
servation during the financial crisis is that one default event may have impact on other 
remaining firms and often causes important losses on the asset values of its counterparties. 
We shall propose a contagion model to take into consideration this phenomenon. 



4.1 Preliminaries 



We begin by generalizing Lemmas 12.11 and 12.21 to the case of processes. Denote by Of 
(resp. Vf) the optional (resp. predictable) cr-field on £1 x K + . 
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Lemma 4.1 1) Any G-optional process Y can be written as Yt = TL A iY t (tj) where 
Y (•) is an Of <S> B(J$J + ) -measurable function on £1 x R + x 

£j ^4m/ G -predictable process Y can be written as Y t = Y^i ^-a 1 ( t i) where Y 1 ^) is a 
Vw <S> £>(IR+) -measurable function onflx M.+ x IR^_. 

Proof. 1) It suffices to consider Y = Zl| s oo |, Z being a ^-measurable random variable. 
By Lemma [2. H for any Jc9, there exists an T s <g> ^(]R^_)-measurable function Z 1 ^) such 
that 

Z = ]T IaiZ^ti). 

ice 

We define 

y/( s a . = / Sjc/ zJ ( s M[a, s y{sj) Uiel\J !{*<**<*} if * > s, 
* /J ' \o if i < s. 

Notice that the process Y^si) is right continuous for any sj. Hence Y 1 is Op <8> 0(Rj_)- 
measurable. By the equality 

II 1 {s<n<t})^[0,s]j( T j) = ^AJ 
IdJ iel\J 

which holds for any J C Q and any t £ [s, +oo[, one can verify that Y t = ]T 7 % a iY i {t i ). 
2) By using Lemma 12.21 a variant of the above argument leads to the predictable version 
of 1). 

□ 



4.2 The model setup 

We consider a portfolio of N assets, whose value process is denoted by S which is an 
^-dimensional G-adapted process. The process S has the following decomposed form 

St = z Z 1 MS I t (r I ), 
ice 

where S[(sj) is Of <8> i?(M+)-measurable and takes values in representing the asset 
values given the past default events tj = sj. Note that St depends on the value of S[ (•) 
only on the set A\, that is, only when t > maxs/. Hence we may assume in convention 
that Sl(si) = for t < s v / where Syj := maxs/ with S0 = 0. 

We suppose that the dynamics of S 1 is given by 
(16) dS{( Sl ) = S[( SI ) * (i4( ai )dt + ^l(8j)dW t ), t > s VI 
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where W is a ^-dimensional Brownian motion with respect to the filtration F, the coef- 
ficients Ht(si) and (sj) are Of ® ^(M^_)-measurable and bounded. Note that for two 
vectors x = (x±, ■ ■ ■ , xjv) and y = (y%, ■ ■ ■ ,vn) in ^ N , the expression x * y denotes the 
vector (xij/i, ■ • ■ ,x N y N ). 

We also suppose that one default event induces a jump on each remaining asset in the 
portfolio. More precisely, for any 1^0, let 

(17) Si vI ( Sl ) = S£,_(*j) * (1 - tS*(*j)), 

where A; = min{i € J|sj = s v /}, J = I \ {k} and j J,k (-) is "Pf <8> /3(M+)-measurable, 
representing the jump at Tyi = Syi given the past default events tj = sj with the last 
arriving default r^. The specification of and J is for treating the case where several 
components in sj are equal to syj. We remind that this convention is harmless since we 
have assumed that Tj 7^ Tj a.s. if i ^ j. 

The value of S-f (•) is determined by Sj. (■) and the coefficients /U^(-) and From 
a recursive point of view, S^(-) actually depends on the initial value So and the coefficients 
indexed by all J C I, together with the jumps at defaults. 

Remark 4.2 We do not suppose that F is generated by the Brownian motion W, allowing 
for some further generalizations. Furthermore, we do not specify the set of assets and the 
set of defaultable names, which permits to include a large family of models. 

We give below several examples. 

Example 4.3 (Exogenous portfolio) We consider an exogenous investment portfolio, that 
is, the underlying assets in the investment portfolio are not included in the defaultable 
portfolio. In this case, the default family contains often highly risky names while the 
investors prefer to choose assets in relatively better situations. However, these assets are 
influenced by the defaults. This is the case considered in [TT] . 

Example 4.4 (Multilateral counterparty risks) The defaults family and the assets family 
coincide, each underlying name subjected to the default risk of itself and to the counter- 
party default risks of the other names of the portfolio. For each name i £ 0, denote by 
S l its value process and by Tj its default time. We suppose that the value of 5* drops to 
zero at the default time Tj, and at the default times Tj where j 7^ i,j € 0, the value of S l 
has a jump. So S z has the decomposed form 

where S^' (•) is J~t <S> S(M^_)-measurable. For I C such that i $ I, let the dynamics of 
S 1 ' 1 satisfy 

dsY(si) = sYis^/is^dt + oY{ SI )dwi), t > SvJ 



13 



where (W 1 ,--- ,W n ) is an n-dimensional Brownian motion with covariance matrice £, 
nl' 1 '(•) and al' 1 (•) are Of ® B{R I + ) -measurable and bounded. In addition, we suppose that 
for any 1/0, 

5^( S/ ) = s^_M(i-7^ fe M) 

where = min{i £ I|sj = Syj}, J = I \ {k}. 

4.3 A recursive optimization methodology 

In this subsection, we consider the optimal investment problem in the contagion risk model. 
Following the recursive point of view on the successive defaults in [5], we have proposed in 
|llj a two-steps — before-default and after-default — optimization procedure. By using 
the density approach introduced in [5], the optimizations are effectuated with respect 
to the default-free filtration F instead of the global one G. This methodology has been 
adopted recently in [13] to add a random mark at each default time by using the joint 
density of the ordered default times and the marks. We discuss this case as an application 
of constrained optimization problems (Remark I4.1U|) using auxiliary nitrations. 

We apply now this recursive optimization methodology to the contagion model de- 
scribed previously. We shall assume the density hypothesis 13.31 in the sequel with v T , the 
law of r, being the Lebesgue measure. Namely, we assume that the random measure of r 
is absolutely continuous with respect to P ® ds on u-field Tt <8> B(W± ), t € M+, where ds 
is the probability law of r. Denote by at(s) the density of \x r on Tt ® B(W±). By taking 
a suitable version, one has that a(s) is an F-martingale for any s 6 R" . Note that we do 
not suppose that the defaults are ordered. 

Let us consider an investor who holds a portfolio of assets, each one subjected to 
contagion risks. The value process of the assets is supposed to satisfy (fl~6|) and (fT7|). The 
wealth process of this investor is described by a positive G-adapted process X and the 
allocation of the portfolio is chosen by the criterion of maximizing the utility of the terminal 
wealth at a finite horizon T. So we are interested in the optimization problem E[{7(Xr)] 
where U is a utility function satisfying the conditions: strictly increasing, strictly concave, 
C 1 on M + and the Inada conditions f/'(0+) = oo, U'(oo) = 0. 

The portfolio of assets is characterized by a G-predictable process tt, representing the 
proportion of the wealth invested on each asset. By Lemma 14.11 the process tt has the 
decomposed form 7r t = J2i ^-a 1 n t ( T i) where (•) is T\®B(M+)-measwc&ble, representing 
the investment strategy given the defaults tj. Hence, to determine an investment strategy 
7r is equivalent to find a family {n 1 (-))icQ- 

Similarly, the wealth process X has the decomposed form Xt = ^2 1 a iX^(ti) where 
X T (-) is Of <8) 0(R+)-measurable, representing the wealth given tj. In view of (fl~6|) . the 
wealth process dynamics is given by 

(18) dXl( Sl ) = Xl{ ai )4{ 8l ) ■ {^ t {si)dt + Ztis^dWt), t > s VI . 
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Note that for two vectors x = {x\, • • • , x^) and y = (j/i, • • • , un) in Mr, the expression 
x • y denotes x\y± + • • • + x^yN- By (fT7|) . we have for any 1^0, 

(19) = X/ yj _( Sj )(l - < 7 ( SJ ) • 7j J5), 

where A; = min{i G I\si = S\ji} and J = I \ {k}. To ensure that the wealth is positive, 
we need to suppose that the cumulative (proportional) losses caused by one default on all 
remaining names is smaller than 1. 

This following simple result is useful for the recursive optimization. 
Lemma 4.5 For any T € M+, 

(20) E[U(X T )] = V / E[U(X^( Sl ))a T (s)]ds. 

f^J[o,T]'x]TM IC 

Proof. We use the decomposed form of X? in Remark 12.81 and take iterated conditional 
expectation to obtain 

E[U(X T )\ = J2n^ T U(xUri))] = Y j ^Wa{U{X i t {t i ))\^ t ]]. 
ice ice 

The lemma then follows by definition of the density. □ 

We introduce the admissible strategy sets. 

Definition 4.6 For I C Q and sj G [0, T] 1 , let A 1 (si) be the set of F-predictable processes 
tt 1 (sj) such that the following two conditions are satisfied: 

!) So Wt{siWt{*i)\ 2 dt < oo; 

2) in the case where 1^0, for any i ^ I and any Si £]s\/i,T], one has 7rf (sj) ■ j£f < 1. 

Denote by A = {(n 1 (-))ice} the set of strategy families tt = (tt 1 (-))ice, where tt 1 (■) 
is a "Pf <8> S(]R^_)-measurable function such that for any si G [0,T] 7 , the process 7r J (s/) is 
in A^si). We say that tt is admissible if tt = (tt 1 (-))j c q g A. 

For our recursive methodology, it will also be useful to consider all the strategies after 
the defaults tj = sj. We introduce the corresponding admissible sets below. 

Definition 4.7 For any / C and any si G [0,T] 7 , let A Dl (si) be the set of families 
(tt k (si, •)) kdi, where tt k (si, ■) is a 7 7 F®^(^+^)-measurable function such that tt k (sk) £ 
A K (s K ) for any s KV G [s v /,T] A V. 

We define the set A^^si) in a similar way. Note that any family tt d1 in A Dl (si) can 
be written as (tt 1 ,^ 1 ) with tt^ 1 G ^2 / ( Sj ) and tt 1 G A 1 (si). 
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Let us now consider the maximization of the utility function on the terminal value of 
the wealth process 



(21) 



J(x,ir) :=K[U(X t )] Xo=x . 



We shall treat the optimization problem in a backward and recursive way. To this end, 
we introduce some notations. Let 



J e (x, s, 7T @ ) := E[U(X®(s))a T (s) | F Sve ] x e 



and 



Vq(x,s)= esssup Jq(x,s,it®) 



We define recursively for I C G, 



Ji{x, si, ir 1 ) := E 



(22) 



U{X^(si)) [ a T {s)dsic 

J]T,+ O o\ ic 



+ E 



}svi,T] 



T.+ocf 

v iu{i} (X^ {t} (siu{i} ) i siu{i}) dsi 



S\/I 



and correspondingly 
(23) 



Vi(x,si):= esssup Ji(x, si, tt 1 ). 



Remark 4.8 Note that viewed from the initial time t = 0, the value of Xl,(sj) depends 
on Xq and all strategies tt j , J C I. However, viewed from the last arriving default r v /, 
its value depends only on the strategy ir 1 if the value of (sj) is given. 

The above constructions provide us a family of optimization problems (Vi(x, s/))/ C 0- 
Notice that at each step, the optimization problem Vi involves the resolution of other 
ones Vj u n\. The whole system need to be dealt with in a recursive manner backwardly, 
each problem concerning the filtration F and the time interval [s\ji,T]. By resolving 
recursively the problems, we can obtain a family of optimal strategies (7? (-))icO- The 
following theorem shows that the global optimization problem, which consists of finding 
the optimal strategy tt £ A for (|2ip is equivalent to finding (vf/(-))i"c8- 

With the above notations, one has in particular 



J (x,7r ) = E[[/p4) [ a T (s)ds + y^f V {i} {X^[ai),8i)da. 

1 J]T,co[ n ■ , J]0,T] 



X -- 



and Vq(x) = sup Jg(x,7r 
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Theorem 4.9 Suppose that V (x, sj) < co a.s. for any I C @, any x > and si £ [0, T] 1 , 
then 

(24) sup J(x,vr) = V 9 (x). 

7TG-4 

Proof. For any / C 0, we deduce from a backward point of view and introduce 

(25) J I (y,s I ,7r Dl )=E\ V / / C/(Xf( S ^))a T ( S )& /c |^ v/ 



KDI 



where Xl vI (sj) = y and 7r D/ is an element in A Dl (sj) (see Definition I4.7p . Note that the 
value of Jj(y,sj,TT Dl ) depends on y, sj and on the choice of ir K (tK) with K D I and 
tj = sj. We shall prove by induction the equality 

(26) esssup Jjiy^^n 31 ) =V I (y,s I ). 

By Lemma H3J we have J^{y^ir) = E,[U(XT)]x =y So the particular case of (|2"6"|) when 
I = is just what need to be proved. 

We proceed by induction on / and begin by I = 0. Observe that Je(y, s, 7r D0 ) = 
Je(y, s, 7r e ). Hence the equality (|2"6"1) holds true by definition when I = Q. Let I be a 
proper subset of O. Assume that we have proved (|26p for all K ~D I. We claim and show 
below that 

J 7 (y, S/ ,7r 3/ ) =E[f7(x£( S/ )) / a T (s)d Sl c 

(27) 7lT '°° [/C 

+ E/ ^um^^/um^^^ r , • 

In fact, by (|25p . the second term in the right-hand side of (|27p equals 
^ / ^E[J /uW (X^«( S7uW ), S/uW ,7r 3/u «)|^ VI ] 



E E / / / E[C/(Xf(^))a T ( S )|^ Sv/ ]d S/ c 

EE/ / / E[C/(Xf(^))a T ( S )|^ VI ] C Z S/ c. 



KDIi£K\I 



We consider, for any i E K \ I, the set ]s v /, T]W x]sj, Tp^W). Note that the subsets 
of Jsv/)^ 1 ]^^ 7 °f the following form 

:= {%\/ I V j G if \ ( J U {<}), s v / < ^ < < T}. 
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are disjoint, in addition, the set 

is negligible for the Lebesgue measure. Hence we obtain 

£ / ds t E [J IU{1} (X 1 ^} ( S/u{4} ) , Slu{t} , ) |^ svj ] 

jg/c J]s V i,T\ 

= Y J I I E[U(X«(s K ))a T (s)\F SvI }d Sl c, 

^jJ]T,oc[ KC J]^i,T] K ^ 

and hence (|27p is established. 

By the induction hypothesis, one has 

(28) esssup Jiu{i}(y,s IU {i},n Dlu{l} ) = V I(J{i} (y,s IU{i} ). 

T D/U{.} e ^/U {i }( Sju{i}) 



Hence we have by (|2< 

esssup J/(y,s/,vr D/ ) < / E[C/(X^(s/))a T (s)|jr Sv/ ]ds/c 

7r3^e^/( S/ ) J]T,oo[-r c 

+ 2 /" E[F 7uW (^«( Sju{i} ) )Sju{i} )|^ Sv/ ]^, 

which, together with the definitions (|22p and (|23p . implies 

esssup Jj(y, s/, 7r D/ ) < V/(y, sj). 

We still suppose the induction hypothesis (f28|) for the converse. For i £ I, s/u{?} G 
[0,T] /LJ W, e > and z G R, there exists a family tPM* } G *A D/uW ( s /u{i}) such that 



/uOl^.s/ui!}.^ 1 ) > ^/u{i}( z ' s /u{i}) 

We fix s J , Xg (sj) = y and 7r J G A^si). By a measurable selection result, we can choose 
a family tt^ 1 G such that 

where i, Zi and s K \j satisfy the following conditions: 

(1) i G I c , s\/j < Si and s, < Sj for all j £ K\(I U {i}), 

(2) ^ = X^«( SjuW ). 
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=)/u W) only depends on z%, Sju{i} and the choice 



— e. 



This implies, by comparing (|22|) and (|27|) . that 



Jj(y, s/, (tt 7 , tt2 j )) > J 7 (y, sj, vr 7 ) - eT|/ c 



By taking the essential supremum over ir 1 £ „4 7 



(sj), one obtains 



7T 



esssup J/(y, s/, tt d ) > V/(y, sj) - eT\I c 



Since e is arbitrary, we get the inequality 



esssup J/(y,s/,7r 3 ) > 



We hence established the equality ([26]) . 



□ 



The existence and the explicit resolution of the optimization problems (Vi(x, si))i c b 
will be discussed in detail in a companion paper. 

Remark 4.10 It is often useful to consider strategies with constraints. In this case, the 
admissible set A (si) for 7r 7 (s/), si £ [0, T] 1 is a subset of A^si) and the admissible 
trading strategy set A for tt is defined similarly : A = {{n 1 ('))ic@} such that for any 
si € [0,T] 7 , 7r / (s/) G A {si). We can also define A^ I (si) and A^ I {si) in a similar 
way. Note that Theorem 14.91 still holds for the constrained strategy. More precisely, let 
us introduce in a backward and recursive way Jj (x, sj, tt 1 ) similarly as in ([22]) and let 
Vf(x,si) = esssup^/g^o Jj(x,si,n I ). Then 



As an application, we consider the case where the losses at default times r = (t\, ■ ■ ■ , r n ) 
are associated to some ^-measurable random variables L = (L\, • • • , L n ). This is the case 
studied in [13] supposing r is a family of ordered default times and using the joint den- 
sity of (r, L) with respect to the default-free filtration. We recall briefly this model and 
show that it can be considered as a constrained optimal problem mentioned in the above 
remark. 

Example 4.11 Let F° be a default-free filtration. The default information contains the 
knowledge on default times n, (i £ 6), together with an associated mark Li taking values 
in some Polish space E. So the global market information is described by the filtration 



sup J(x,tt) = V£(x). 



G° = F° V D 1 V • • • V B" V B Ll V • • • V D Ln , 
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where for any i 6 6, the filtration B* is as in Section \2. 11 and the filtration D^* = {T>^ l )t>Q 
is defined by T>^ 1 = cr(Ljl{ r .< s }, s < t) made right-continuous. Note that any G°-optional 
(resp. predictable) process Z can be written in the form 



z t = Yl 1 A I t Z t('>~i, L T ), (resp. Z t = \ a i_Z\{t u L/), ) 




where Z 1 '(•) is a F ° ® <S(R+ x M J ) (resp. P F ° <8> B(K+ x JR 7 ))-measurable function, Lj = 
(Lj)j g /. In particular, the control process ir can be written as 



We explain below how to interpret the above model as a constrained optimization prob- 
lem. The point is to introduce suitable auxiliary nitrations. Let F := F° V o~(L\, ■ ■ ■ , L n ). 
It is the initial enlargement of F° by including the family of marks L. Define also 
F 7 := F° V a(Li,i £ I). Observe that the control vr 7 (-,L/) is actually V v i ® B(R\)- 
measurable and is hence "Pf <8> S(]R^_)-measurable. The wealth X 1 (•, Lj) is O w i <S> £>(M+ )- 
and hence Of <g> i3(M / f )-measurable. We introduce the filtration G = F V D 1 V • • • D n , which 
is the progressive enlargement of the filtration F with respect to the default nitrations. 
Note that G D G°. 

The Example 14.111 can be considered as a constrained problem by using the auxiliary 
nitrations F and G. Indeed, an admissible control process n has now the decomposed 
form ir t = J2ic0 ^A 1 7r t( r /)) where 7r 7 (s/) is "Pf ® S(K+ )-measurable and 7r / (s/) € ^4 7 (s/), 
the admissible set A (sj) being defined in Definition 14.61 Let us now make precise the 
constrained admissible strategy sets: let Ai(si) be the subset of A 1 (sj) such that ir 1 (sj) 
is Vfi <8> S(M^_)-measurable if 7r J (s/) G ^(sj). By Remark I4.1UI we can apply Theorem 
14.91 to solve the problem. We finally remark that we only need the density hypothesis of 
r with respect to the filtration F but not necessarily the stronger one on the existence of 
the joint density of (r, L) with respect to F°. 
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